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SEARCHING FOR NON GAUSSIAN SIGNALS IN THE BOOMERANG 2003 CMB MAPS 
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ABSTRACT 

We analyze the BOOMERanG 2003 (B03) 145 GHz temperature map to constrain the amplitude 
of a non Gaussian, primordial contribution to CMB fluctuations. We perform a pixel space analysis 
restricted to a portion of the map chosen in view of high sensitivity, very low foreground contamina- 
tion and tight control of systematic effects. We set up an estimator based on the three Minkowski 
functionals which relies on high quality simulated data, including non Gaussian CMB maps. We find 
good agreement with the Gaussian hypothesis and derive the first limits based on BOOMERanG data 
for the non linear coupling parameter /nl as —300 < /nl < 650 at 68% CL and —800 < /nl < 1050 
at 95% CL. 

Subject headings: cosmology: cosmic microwave background 



1. INTRODUCTION 

While cosmology is entering its precision era, the tar- 
get of experiments aimed at the Cosmic Microwave Back- 
ground (CMB) is shifting towards weak signals, e.g. po- 
larization, the SZ effect and non Gaussian (NG) features. 
Detection of NG signals can be of significant help in 
constraining the mechanisms that explain the genera- 
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tion of cosmological perturbations. Provided that sys- 
tematic effects will not degrade high sensitivity CMB 
mapping, present and future experiments could in prin- 
ciple be sensitive to non linearitie s due to second or - 
der effects in perturbation theory (jBartolo et al. 20041 ) . 
This signal is usually parametrized by a non linear cou- 
pling factor /nl that controls the level of a quadratic 
contribution to the primord ial gravitational potential $ 
(|Komatsu and Spergel 20011) : 

$(x) = $ G (x) + /nl [<Mx) 2 - (<Mx) 2 )] (1) 
where $g is a zero mean, Gaussian random field. 

Several groups have reported NG constraints on 
CMB data. All suborbital efforts to date have 
found no significant deviation from Gaussianity in 
the C MB field: MAXIMA- 1 reported |/ NL < 950 
at la (|Santos et al. 20031 : IWu et al 200H), while VSA 
found an upper limit of 5400 at 2a (|Smith et al. 20041) : 
Archeops has recently improved their limits to —800 < 
/nl < 1100 (2a), although their analysis is based on as- 
sumptions only valid for the large angular scales dom - 
inated by the Sachs- Wolfe effect (jCurto et al. 2 006). 
The BOOMERanG 1998 data set has also been t ested 
for Gaussi anity, both in pixel ( Polenta et al. 20021 ) and 
in Fourier (|De Troia et al. 20031 ) space, finding no trace 
of NG signals. However, BOOMERanG has set no /nl 
limit so far. One of the purposes of this paper is to pro- 
vide such limits with the analysis of the new 2003 data. 
The limits presented here are more stringent than those 
found by any suborbital experiment to date, properly 
accounting for sub-horizon angular scales. 

T he WMAP team con strained /nl to be -54 < / N l < 
114 ()Spergel et al. 20061 ). Using an improved version of 
the WMAP team estimator Creminelli et al. (2007) set 
the most stringent limits to date at —36 < /nl < 100. 
Thus the /nl analysis does not show any departure from 
Gaussianity in WMAP data. However, some authors 
have looked at general deviations from Gaussianity (i.e. 
not based on any specific parametrization of NG) and 
claimed highly-significant detection of NG features in 
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the WMAP map s (jCopi et al. 2004 IVielva et al. 2004 
iCruz et al. 2001 . 

In this paper we perform a pixel space analysis of 
the B03 temperature (T) data set, using the observed 
field's moments and Minkowski functionals (MFs) to 
build Gaussianity tests. We assess the statistical sig- 
nificance of our results comparing the data to a set of 
highly realistic, Gaussian Monte Carlo (MC) simulated 
maps. In order to constrain /nlj we build a goodness of 
fit statistics based on MFs and calibrated against a set 
of NG CMB maps, th at are generated ac cording to the 
algorithm set forth in (Lig uori et al. 2003f ). 

The plan of this letter is as follows: in section [2] we 
briefly describe the B03 experiment, the dataset it has 
produced and our simulation pipeline. In section [3] we 
compute the map's moments and MFs of the data and 
compare results against Gaussian MC simulated maps. 
Then we derive constraints for /nl- Finally, in section [4] 
we draw our main conclusions. 

2. THE BOOMERANG 2003 DATASET 

The balloon borne B03 experiment has been flown from 
Antarctica in 2003. It gathered data for 14 days in three 
frequency bands, centered at 145, 245 and 345 GHz. Here 
we restrict ourselves to the 145 GHz data that are most 
sensitive to CMB fluctuations. These have been obtained 
with polarization sensitive bolometers (PSB). The anal- 
ysis of the dataset has produced high quality maps of 
the southern sky that have been conveniently divided in 
three regions: a "deep" (in terms of integration time) 
survey region (~ 90 square degrees) and a "shallow" sur- 
vey region (~ 750 square degrees), both at high Galactic 
latitudes, as well as a region of ~ 300 square degrees 
across the Galactic plane. The deep region is completely 
embedded in the shallow region. Here we only consider 
a subset of the data that contains all of the deep region 
and part of the shallow, for a total of 693 square de- 
grees (1.7% of the sky). The mask we use is square, 
26 degrees in side, centered at about RA=82.6° and 
DEC=— 44.2°, and excludes all detected point sources in 
the field. This region has been selected in view of high 
sensitivity CMB observation with low foreground con- 
tamination and was observed with a highly connected 
scanning strategy to keep systematics under control. We 
use t he T data map red uced jointly from eight PSB at 145 
GHz (Masi et al. 2005T ). In this region, the signal rms on 
3.4' pixels is ~ 90 /iK and instrumental noise has an rms 
of ~ 20 /iK in the deep region and of ~ 90 /iK in the shal- 
low region. In harmonic space, binned estimates of the 
CMB angular power spectrum retain signal to noise > 1 
well beyond I ~ 1000. One may compare these figures 
with WMAP: in the three year release, WMAP combined 
sensitivity in the region observed by B03 is ~ 100 /iK on 
3.4' pixels, close to WMAP's mean pixel error. However, 
WMAP's beams are broader than B03, so its I space er- 
ror is ~ 5 times larger than B03 at t ~ 1000. On the 
other hand, B03 has not been devised to measure mul- 
tipoles at I < 50. In this sense, our NG analysis probes 
angular scales complementary to those constrained by 
WMAP. While we do not consider here the Stokes Q 
and U polarization maps, our T map has been marginal- 
ized with respect to linear polarization. For a descrip- 
tion of the instru ment and the mea sured T and polar- 
ization maps see (|Masi et al. 20051) and for the CMB 




Fig. 1. — The left panel shows the distribution of S3 estimator 
calculated from the 200 Gaussian MC maps. The S3 value of the 
B03 data is represented by the vertical line. The right panel shows 
the same for S4 

TT, TE and EE power spectra see (iJones et al. 20061 
iPiacentini et al. 20061 : iMontrov et al. 2006D . 

To assess the robustness of our tests of Gaussianity 
we use a set of simulated MC maps that mimic the B03 
data. To produce these, we follow the same steps per- 
formed when analysing real data. The Gaussian CMB 
sky signal is simula ted from the power spe ctra that best 
fits the B03 data (jMacTavish et al. 2006h . This signal 
is smoothed according to the measured beam and syn- 
thetiz ed into a pixelized sky map, using Healpix rou- 
tines (jGorski et al. 2005V Using the B03 scanning strat- 
egy, the signal map is projected onto 8 timestreams, 
one for each 145 GHz detector. Noise only timestreams 
are also produced, as Gaussian realizations of each de- 
tector's noise power spectral density, estimated from 
the data accounting for cross talks among detectors. 
The t imelines are reduc e d with the ROMA mapm aking 
code (jNatoli et al. 20011 : iDe Gasperis et al. 2005? ) repli- 
cating the actual flight pointing and transient flagging, to 
produce T,Q and U maps. With this procedure, we can 
simulate signal, noise and signal plus noise timestream. 

To constrain /nl we use MC simulations of NG CMB 
maps obtained from a primordial gravitational potential 
of the form given in eq. ([lj. These maps have been 
produc ed including first o rder CMB radiative transfer 
effects ([Liguori et al. 20031) . The power spectrum of the 
NG maps is identical to that of the Gaussian CMB sim- 
ulations. 

3. TESTS OF GAUSSIANITY AND CONSTRAINTS ON / NL 

Working at 6.8' Healpix resolution (A^idc = 512), we 
first compute the normalised skewness S3 and kurtosis 
S4 of our pixelized field Tj. These are obtained from the 
variance a 2 = 1/{N - l)J2ii T i - ( T )) 2 and from thc 
third and the fourth moment /i3 = Ylii^i ~ (T)) 3 /N an d 
= T,i( T i - (T))*/ 1 ^, wh ere TV is the total number 
of pixels of the map and (T) — ~^2nTi/N its mean. We 
have S3 = /i3/cr 3 , S4 = /i4/fT 4 — 3. From the data we 
get 53 = —0.016 and S4 — 0.096. These values are plot- 
ted in Fig. [T] as a vertical line and compared to the em- 
pirical distribution as derived from the MC (signal plus 
noise) maps. From the latter we compute the probabil- 
ity P{Sf m > S , 3 Jata ) = 58% and P(Sf m > Sf ata ) = 77%. 
Hence, for these tests the data are compatible with the 
Gaussian hypothesis. The same tests are repeated after 
having degraded the map to 13.6', finding similar res ults. 

To analyse the map with MFs (|Gott et al. 19901 ) we 
consider the excursion sets Q defined as the map's sub- 
sets exceeding a given threshold v. Q(v) = {Ti : (Tj — 
(T))/a > v}. The three MFs measure the total area of 
the surviving regions of the map (-Mo), their total con- 
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Fig. 2. — Left panels show the MFs estimated from the B03 data 
(diamonds) and the 2<r confidence limits from the 200 Gaussian MC 
maps. From top to bottom, the results correspond to the MFs MO, 
Ml and M2, respectively. Right panels show the residuals between 
the MC mean and the data. 

tour length (Mi), and the genus of the distribution which 
is related to the difference between the number of "hot" 
and "cold" regions (M 2 ). For a Gaussian field the ex- 
pectation values of the functionals depend on a single 
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In the case 



parameter r: (M (i/)) = | 1 — erf 

#exp(-if), (M 2 (v)) = ^exp 
of a pure CMB signal (no noise), r is given by r 
j Y,^ (21 + 1) 1(1 + (ISchmalzing and Gorski 19911 
iWinitzki and Kosowskv 19981 ). Hence, M x and M 2 de- 
pend on the power spectrum Ci. It is hence critical 
that the simulations reproduce the model CV's that best 
fits the data. We work in the flat sky limit, projecting 
our T values on th e plane locally tangent to the map 
(jCabella et al. 2004f) . In Fig. [2] we plot MFs for the data 
and 2a limits set by 200 Gaussian simulations, as well as 
the data residuals and their (again, 2a) errors. 

The functionals are computed at 9 thresholds evenly 
spaced between —4a and +4er. 

Using the MC maps we can define a x 2 statistic: 



Xb, 



(Mf - {Mf m )) v CrL(M? - (Mr 1 )),'- (2) 



Here Mf (Mf m ) is any of the three MFs obtained from 
the data (simulations), (•) is the mean value over MC re- 
alizations, andC^,,/ = ((M iiV - (M itV ))(Miy - (Miy))) 
is a covariance matrix, estimated from an independent 
set of ~ 1000 Gaussian maps. In the first three panels 
of Fig. [3] we show xf for each MF (vertical line), plotted 
along with the empirical distribution sampled via MC. 
We can define a "joint" estimator by grouping the Mi's 
in a single, 27 elements data vector Mj = {Mo, Mi, M 2 } 
and defining a Xb as a trivial extension of Eq. [2] It is 
important that the covariance matrix built for the joint 
estimator correctly accounts for correlations among dif- 
ferent functionals. However, we have found that in order 
to pin down to a stable level these off-block couplings, one 
requires a number realizations significantly higher than 



Fig. 3. — The \% distribution (histogram) of MC simulated maps 
and data value (vertical line) for the MFs. Top: area and contour 
length. Bottom: genus and "joint" estimator (see text). 



the ~ 1000 used throughout our analysis. The latter 
number cannot be realistically increased to desired level, 
because the GLS map making procedure is a demand- 
ing computational task, even for the supercomputers we 
have used. Fortunately, we have found that using white 
(instead of correlated) noise to estimate the covariance 
matrix has a negligible impact on the analysis. This find- 
ing can indeed be justified a posteriori, because the GLS 
map making procedure is very effective in suppressing 
noise correlations, that contribute very weakly to the es- 
timator's final covariance. The joint x 2 of the data is 
displayed as the fourth panel in Fig. [3j along with the 
MC empirical distibution. The probability P(x 2 > Xb) 
that a Gaussian map has a larger x 2 than the B03 map 
is 76% for M , 83% for Mi, 76% for M 2 and 67% for the 
"joint" estimator. The values are fully consistent with 
the Gaussian hypothesis. We conclude that our pixel 
space analysis does not detect any sign of NG behavior 
in the B03 data. 

We now want to constrain the quadratic coupling pa- 
rameter /nl defined in Eq. [TJ To this purpose we simu- 
late NG CMB realizations in the following way: firstly, 
we generate the Gaussian and NG part of the primor- 
dial potential defined by Eq. ([1]) ; then we convolve them 
with CMB first order radiation transfer functions to get 
the final CMB sky. In this way we produce 200 G 
(Gaussian) maps and 200 NG counterparts (each G map 
has a uniquely defined NG counterpart), so that for a 
given /nl our (signal only) map is G + /nl * NG. By 
adding noise maps, we can define MF estimators in the 
spirit of the section [3l with the difference that they are 
now functions of / NL : J b (/nl) = Mf - (Mj(/ NL )) 
(we only consider the "joint" estimator in what fol- 
lows). Consequently, we now define the data x 2 as 
Xb(/nl) = Jb(/nl) t C _1 J b (/ N l)- While in principle 
the covariance of the Mj's is a function of /nl, this 
depende nce is expected to b e weak and is usually ne- 
glected (|Komatsu et al. 20031) . We have tested for this 
explicitely by using our NG simulations. We plot x% 
as a function of /nl m Fig. |4] (left panel). Goodness 
of fit analysis yields -300 < /nl < 650 at 68% CL 
and -800 < /nl < 1050 at 95%, with a best fit value 



(x B at its minimum) of /nl = 200. In order to cross 
check this result, we estimate a "frequentist" confidence 
interval for /nl by sampling the empirical distribution 
of the Mj-based x 2 -, computed for / NL = 200. The 
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Fig. 4. — The left panel shows the B03 data x 2 (x% m tne 
text) as a function of /nl, while the panel on the right shows the 
empirical x 2 distribution for NG maps (/nl = 200), sampled using 
1000 simulations 

resulting histogram is displayed in the right panel of 
FiglH This analysis shows that we should expect to con- 
strain |/nl| ;$ 1000 at 95%, thus suggesting that our 
limits derived through goodness of fit analysis are con- 
sistent. One may consider what limits on /nl would be 
derived if we use, in place of MFs, the map's skweness 
and kurtosis defined in Sect. [3] as elements of a two di- 
mensional data vector. We thus repeated our goodness of 
fit analysis using these statistics and found weaker limits: 
-950 < /nl < 1350 at 68% CL (-1400 < / NL < 1800 
at 95% CL). Even so, it is quite remarkable that a crude 
1-point pixel statistic degrades the final sensitivity only 
by a factor ~ 2. Of course, in order to find "optimal" 
constraints one has to implement a more sophisticated 
analysis. 

The B03 dataset is not expected to be noise dominated 
at the 6.8' resolution employed in the analysis above. To 



show that this is the case, we repeated all of our proce- 
dures for 13.6', finding similar (though slightly weaker) 
constraints. 

4. CONCLUSIONS 

We analysed the B03 145 GHz T map in search of 
NG signals. We worked in pixel space at 6.8' and 13.6' 
Healpix resolution. We computed the skewness and kur- 
tosis of the map, as well as its three MFs. We compared 
these estimates against a set of simulated Gaussian maps 
that have been reduced using the same analysis pipeline 
as real data, finding no evidence of NG behavior. To 
quantify the latter statement, we define goodness of fit 
statistics jointly based on all three MFs, showing that the 
probability for a Gaussian simulation to have a x 2 larger 
than the data is ~ 67%. Assuming a model for primor- 
dial fluctuations that predicts a quadratic perturbation 
to the gravitational potential, we set limits on the non 
linear coupling parameter /nl as —800 < /nl < 1050 at 
95% CL (-300 < / NL < 650 at 68% CL). These limits 
may be regarded as complementary to the constraints set 
by WMAP in view of the better signal to noise ratio at 
high resolution in the B03 field. 
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